Gravitational and electromagnetic outputs from binary neutron star mergers 
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The late stage of an inspiraling neutron star binary gives rise to a strong emission of gravitational 
waves due to its highly dynamic, strong gravity. Interactions between the stellar magnetospheres, 
driven by the extreme dynamics, can produce considerable outflows. We study the gravitational 
and electromagnetic waves produced during the inspiral and merger of a binary neutron star system 
using a full general relativistic, resistive MHD evolution code. We show that the interaction between 
the stellar magnetospheres extracts kinetic energy from the system and powers radiative Foynting 
flux and heat dissipation. These features depend strongly on the configuration of the initial stellar 
magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries 
and has a distinctive angular and time-dependent pattern. 



PACS numbers: 

Introduction: Binary systems involving neutron stars 
are among the most likely sources of detectable gravita- 
tional waves, making them important sources for gravi- 
tational wave detectors such as Advanced LIGO. Grav- 
itational wave observations of neutron star binaries are 
expected in the next few years, and the combination of 
these data with electromagnetic observations will pro- 
vide new opportunities to study the fundamental physics 
of these stars, and their interactions with their envi- 
ronments. As an example of the important interplay 
between gravitational wave (GW) and electromagnetic 
(EM) data, consider models of compact object mergers 
for short, hard gamma ray bursts (SGRBs) (see e.g. |l] 
for a review). If the initial binary consists of two neu- 
tron stars, the merger may result in either the immediate 
collapse to a black hole or the creation of an intermedi- 
ate, hyper-massive neutron star, supported by thermal 
pressure and differential rotation, followed by a delayed 
collapse to a black hole. The interaction of the central 
compact object with a magnetic field can power radiation 
with the hard spectrum and short time scales character- 
istic of SGRBs. Considerable effort has been devoted to 
understanding these bursts, concentrating on emission 
after the merger. Recently however, there has also been 
a significant interest in possible EM emissions preceding 
the merger. This has been motivated, in part, by both 
the desire to maximize opportunities for gravitational 
and electromagnetic waves detection, as well as surveys 
for EM transients that may be precursors to SGRBs 
(e.g ^) and other transients from compact binary merg- 
ers (e.g. [3]). A basic question here is what possible mech- 
anisms prior to merger can yield sufficiently strong EM 
emissions to be detected. These detection might be aided 



by GW observations, which in future years will provide 
timing information, sky localization, and other physical 
parameters to aid in identifying sources (e.g. [SHS]). 

Precursor emissions may be generated by crust crack- 
ing due to resonance effects [7] or magnetosphere interac- 
tions [8HT2I . A particularly appealing mechanism, unipo- 
lar induction^ relies on the extraction of a star's kinetic 
energy by the magnetically dominated magnetosphere 
surrounding the stars. Here an "effective circuit" is es- 
tablished between the stars, powered by an electromotive 
force produced by the relative motion of a weakly, or non- 
magnetized star through the magnetic field of the other, 
allowing for energy dissipation. The energy released by 
this circuit either reaches its maximum prior to merger 
or gives rise to episodic emissions depending on the re- 
sistance of the circuit. Different models build upon this 
possibility |TQHT2] , predicting a transient preceding the 
merger on a time of a few seconds, possibly accompanied 
by a radio signal. 

The expected luminosity for a binary system composed 
by a magnetized primary star and an unmagnetized sec- 
ondary one is roughly given by [12 , 

Lind ~ 7 X 10^°^{Bt,f (aaokm)"''/' ergs/s, (1) 

where ^ gives the azimuthal twist of flux tubes connecting 
both stars, the magnetic field strength of the primary 
star and we have assumed typical sizes = Rc = 10 km 
and masses = Mc = 1.4 for the stars. A key 
unknown here is the azimuthal twist which depends on 
the total resistance of the circuit. For a free- space mag- 
netosphere ^ = 647r'i;rei, where Vrei the relative veloc- 
ity between the stars, and the resultant luminosity is 
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Lind ^ 1.7 X 10"^^ {Bli)'^a~^^^ ergs/s. The above es- 
timate assumes a quasi- adiabatic process, with an un- 
magnetized companion so that the induced circuit can 
be analyzed in simple terms. However a richer structure 
should arise when both stars are magnetized; in partic- 
ular, significant reconnect ions of the magnetic field lines 
can take place, which would depend on the orientation 
of magnetic moments. Estimates of these are difficult 
to obtain, requiring numerical solutions to unravel the 
possible phenomenology. 

To study EM emission from dynamic neutron star sys- 
tems, it is crucial to properly model global electromag- 
netic fields, as well as gravitational and hydrodynamic ef- 
fects. That is, to consistently incorporate magnetic fields 
present in the stars, their effects in the magnetosphere 
region surrounding the compact objects (see discussion 
in [l3]), the rapidly changing, and extreme, gravitational 
fields involved, as well as the magnetized matter describ- 
ing the stars. We have done so using general relativity 
coupled to relativistic, resistive magnet ohydrodynamics 
in full 3D. Details of our implementation and work build- 
ing up to it can be found in |13H2Q] . 
Approach and techniques: We study inspiraling, 
magnetized binary neutron stars, concentrating on the 
pre-merger epoch, and monitor closely the EM and GW 
radiation produced by the system, the ensuing current 
sheets, dissipation layers, and reconnect ion regions. We 
focus on a binary initially described by a pair of equal 
mass neutron stars in a quasicircular orbit, with an ini- 
tial separation L = 45 km, corresponding to an orbital 
frequency flo = 1-85 rad/ms. The geometric and mat- 
ter initial data for this system are obtained with the 
LORENE library [21], assuming a polytropic equation 
of state P = Kp^ with T = 2 and = 123. The indi- 
vidual stars have a baryonic mass M = 1.62 M©, radius 
= 13.6 km, and a magnetic moment fii that pro- 
duces a dipolar magnetic fields in each star's comov- 
ing frame. The magnetic moment is related to the radial 
magnetic field at the pole of the star by the relation 
jii = B^R^. In our simulations we set 5* = 1.5 x 10^^ G, 
a values on the high end of observations from binaries 
but still realistic. The electric field is obtained from the 
ideal MHD condition E = — v x B, where the velocity in 
the star is given by the orbital motion, and we assume 
the magnetosphere is initially at rest. Our numerical do- 
main extends up to L = 320 km and contains five nested 
fixed mesh refinement (FMR) grids, each finer grid with 
twice the resolution of its parent grid, such that the high- 
est resolution grid has Ax = 0.45 km and extends up to 
58 km, covering both stars and the inner part of the 
magnetosphere. We compute the Poynting- vector lumi- 
nosity at three different surfaces, the furthest located at 
R^^^ 180 km. 

We consider three related initial configurations of the 
magnetic moments to gain insight into the overall be- 
havior of magnetized binaries: aligned and equal mag- 



netic moments (Case U/U) with pi = /j,2 = /i, anti- 
aligned and equal magnetic moments (Case U/D) with 
pi = — /i2 = and aligned magnetic moments with one- 
dominant moment (Case U/u) with jui = 100 /i2 = 
The direction of the magnetic moments is given with re- 
spect to the orbital angular momentum direction. 
Results: In all cases, the magnetic field strengths have a 
negligible effect on the orbital dynamics of the system up 
to merger [HI [22] , as their contribution to the total iner- 
tia is several orders of magnitude below the matter contri- 
bution. Consequently, the three cases we study progress 
to merger in exactly the same way, producing the same 
gravitational signal. The stars orbit each other, and as 
a consequence of gravitational waves carry both bind- 
ing energy and angular momentum, the system merges 
in ^ 2.5 orbits. At this point we stop our simulations 
(leaving the post-merger epoch for future work). 

While the orbiting behavior is the same for all three 
cases considered, the electromagnetic field dynamics and 
magnetosphere interactions depend sensitively on the ori- 
entation of the magnetic dipolar moment of the stars. 
Such interactions affect strongly the resulting field's 
topology and induced dissipation regions, reconnections 
of field lines, and net Poynting flux, as well as other rel- 
evant features. At a basic level, the accelerated orbital 
motion of the stars induces only a small degree of winding 
along its magnetic field lines; thus the magnetospheres es- 
sentially co-rotate with the stars and the magnetic field 
at their surfaces (and therefore, the magnetic dipolar mo- 
ments) remains almost constant until the merger. 

In the anti- aligned case (U/D), a shear layer is in- 
duced at the midplane, separating two regions filled with 
magnetically dominated plasma moving in opposite di- 
rections. The poloidal component of the magnetic fields 
points in opposite directions in the midplane, allowing for 
reconnections that consequently connect both stars (see 
Fig. 1). Further, the binary's rotation stretches these 
magnetic field lines -increasing their tension- and gives 
rise to a toroidal magnetic field that, beyond some dis- 
tance from the orbital center and in the plane connect- 
ing both stars, has equal signs above the equator and the 
opposite sign below. The region near the stars, and espe- 
cially near the orbital plane, is more complicated. Recon- 
nections at the midplane seed a current sheet that prop- 
agates outwards forming a spiral pattern. Interestingly, 
the current sheet structure rotates with the periodicity 
of the orbital motion, effectively providing a "spacetime 
tracer", that would induce distinct features in possible 
EM signals from these systems. 

When the magnetic moments are aligned and equal 
(Case U/U), a shear layer is also induced in the mid- 
plane between the stars due to the relative opposite ve- 
locities of the magnetospheres. (See Fig. 1.) However, in 
this case the symmetry and consequent magnetic repul- 
sion between aligned field lines prevents reconnections 
between the stars. Again, rotation induces a toroidal 
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magnetic field, although the topology is clearly different 
than the one observed in the U/D case. Here the struc- 
ture is reminiscent of the aligned (dipole) rotator [23] . 
This is natural, as there is a net effective dipole at leading 
order with /ieff = 2/i. However in this case the symme- 
try of the system implies an (approximate) periodicity in 
the solution given by half the orbital period. A current 
sheet is also induced at the equatorial plane. This cur- 
rent sheet arises at far distances first (defining the light 
cylinder), and closes in as the orbit tightens, displaying 
a behavior analogous to the current sheet of an aligned 
rotator. Again the binary's dynamics modifies the struc- 
ture and behavior of this current sheet. A spiral pattern 
results from the inhomogeneous structure of the electro- 
magnetic fields produced by the binary structure and the 
orbital motion. 

When the first star has a dominant magnetic field mo- 
ment and the second star is significantly less magnetized 
(Case U/u), the magnetic field of the first star eventually 
dominates that of the companion even near its surface. 
A net magnetic flux from the stronger magnetized star 
threads the weakly magnetized star, and a circuit can be 
established. As a consequence, global electromagnetic 
field features are roughly given by the inspiraling dipole 
together with the induction on the other star. An inter- 
esting effect arises as magnetic fields from the strongly 
magnetized star slide off the companion's surface and re- 
connect afterwards. This produces a dissipation tail as 
illustrated in Fig. 1. This tail gradually grows in extent 
as the merger progresses, populating a current sheet. 

The qualitative understanding of the radiation from 
these three systems, including the angular distribution, 
can be obtained from the Poynting fiux shown in Fig. 2. 
Both the U/D and U/U cases radiate more strongly along 
the shear layer between the two stars, while the U/u case 
does so mainly on the equatorial plane and primarily in 
the direction of the strongly magnetized star. The radi- 
ation in the U/D and U/U cases is partially collimated. 
The flux density in a polar cap (with opening angle of 
9o < 30^) is larger than the average by 2.5 and 1.9, ac- 
counting for 1/3 and 1/4 of the total power, respectively. 
The radiation from the U/u case is emitted mainly near 
the equatorial plane, with 2/3 of the total energy radi- 
ated between 60^ < Qo < 90^. 

The total Poynting luminosity (Fig. 3) in time shows 
several important features. First, the U/D case is sig- 
nificantly more radiative than the U/U case. This is in- 
teresting as the "inner-engine" in both cases is powered 
by the fields of each star and the orbital motion, which 
are the same in both cases (modulo the field topologies), 
implying a more efficient tapping of orbital energy with 
anti-aligned magnetic moments (Case U/D) when they 
are aligned (Case U/U), possibly due to the energy stored 
as toroidal magnetic tension in the U/D case. 

To explain the observed phenomenology, guidance can 
be drawn by considering an effective magnetic moment 




FIG. 1: Magnetic field configurations (field lines) and current 
sheets (orange regions) for the U/D (top panel), U/U (mid 
panel) and U/u (bottom panel) respectively at times t=3.6ms 
and t=4.7 ms. A shear layer is induced in a midplane between 
the stars in the U/D and U/U cases, with reconnections given 
rise to field lines connecting both stars in the U/D (top panel) 
and repulsion between the magnetic fields in the U/U case 
(middle panel). A dissipation tail emanates behind the weakly 
magnetized star (bottom panel). 




FIG. 2: Poynting flux for the U/D (top panel), U/U (mid 
panel) and U/u (bottom panel) case at t = 4.1 ms, in a volume 
(left panels) and evaluated in a sphere located at r = 80 km 
(right panels). 



and how it radiates within a force-free magnet osphere. 
In particular, recall that an aligned (dipolar) rotator 
within force- free radiates a power Ldip ~ jj.'^fl^ [23] . 
Since the U/U case has, to leading order, a rather con- 
stant effective magnetic dipole /ieff = 2/i, one can expect 
L^/^ ^ Q^. This behavior is consistent with our re- 
sults for such case at late times, as can be observed in 
Fig. 3. The anti-aligned configuration U/D is more in- 
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FIG. 4: Currents (drawn as arrows) and charge density (color 
coded) for the U/D (first panel), I 
time (ms) ^ ' (third panel) cases at t — 5.9 ms. 



-^^jq39| 1 , 1 , ± , 1 , ± , ± J-^^-^q53 coded) for the U/D (first panel), U/U (second panel) and U/u 



FIG. 3: The gravitational (right axis) and electromagnetic 
(left axis) luminosities for the three different magnetic field 
configurations, together with two curves with and Q^^^^ 
dependency. 



volved, since the leading order of the effective dipolar 
moment vanishes and higher order terms need to be con- 
sidered. For a rotating quadrupole, Lquad ^ Q^^^, with 
Q a measure of the magnetic quadrupole from the sys- 
tem, which scales as df^ (with d the characteristic size of 
the quadrupolar source). An effective quadrupole could 
be defined here as the stars approach each other and 
fields reconnect. The characteristic size of the effective 
quadrupole necessarily scales with the binary's separa- 
tion (a), though the magnetic field, as sourced by the in- 
dividual stars, remains fairly constant. Thus, one would 
expect Lquad ^ a^O^^ ^ 1^^^^ (using the Keplerian re- 
lation a (X 1^-2/^). This results in a rather flat profile 
within the frequencies spanned by the system in its last 
2.5 orbits, consistent with the observed behavior. 

However, in the U/u case this picture is not quite appli- 
cable. Here the dominant field is given by a non-rotating 
star and, as opposed to the other cases, the contribution 
of both stars does not give rise to an effective rotation to 
define a light cylinder before merger. Nevertheless, the 
unipolar induction model can be adopted. This model, 
using again the Keplerian relation, leads to an estimated 
luminosity L^J^ ~ Q}^^^ . This behavior fits reasonably 
well the observed late time behavior, as shown in Fig. 3. 

Notably, a circuit seems to be established in all cases 
(see Fig. 4), albeit with different characteristics. In both 
the U/U and U/D cases, the circuit is established in be- 
tween the stars and along the axial directions, while in 
the U/u case the circuit closes between the stars. Addi- 
tionally, in the U/u case the amount of energy dissipated 
as Joule Heating {JiE^) is comparable to the Poynting 
energy radiated. For the U/D and U/U cases, on the 
other hand, the energy dissipated as heat is only a frac- 
tion of the radiated energy, between 25-50%. This in- 
dicates that magnetosphere interactions in the U/D and 
U/U cases induce more radiation than can be explained 



via an unipolar induction model. 

Discussion: Our results reveal the late orbiting stages of 
binary neutron stars can induce strong electromagnetic 
emissions, in addition to strong gravitational wave out- 
put. These electromagnetic emissions are sourced by the 
ability of the magnetosphere to tap orbital energy from 
the system and induce both a flux of (Poynting) elec- 
tromagnetic energy as well as current sheets and heat 
dissipation. Indeed, there may be other possible emis- 
sions coming from the current sheets produced in the 
U/U and U/D cases, where strong cooling can give rise 
to gamma-rays [24l [S^l via synchrotron [25] or inverse 
Compton scattering [26 . (See also discussion in [27 ). 
Interestingly, the induced current sheets are heteroge- 
neous and have the orbital behavior strongly imprinted 
on them, which, in turn, could give time-varying signals 
that might aid in detecting these systems electromagnet- 
ically. Furthermore, accelerating fields can arise natu- 
rally at gaps [28H3T] , energizing a population of particles 
that emit high energy curvature and synchrotron radi- 
ation. Just as in the case of pulsar emissions, several 
of the above mentioned mechanisms might be at play 
to yield particular signals. Understanding which mech- 
anism are the most important is an active area of re- 
search. At the simplest level however, if the medium is 
optically thick (as has been argued would be natural for 
the magnetic field strengths considered e.g. [I0l|24]), the 
radiation will produce a relativistic outflow that will emit 
as a black body with an effective temperature given by 
the injection energy, in particular of Poynting luminos- 
ity. Further, recall that magnetic fields will not affect 
the binary dynamics provided the strengths considered 
are < 10^^ G. Thus, our Poynting luminosities can be 
rescaled in strength by (5/10^^ G)^. From our simula- 
tions we can directly extract the total Poynting luminos- 
ity emitted in a radius ^ 30 km (at, e.g. t ^ 5 ms) and 
express the resulting temperature T = Sa/b B\{^ K, 
with the constant Sa/b providing the value for the differ- 
ent cases {Su/u — 4.7, Su = 2 and Su/d = 10)- Notice 
that the luminosity in Case U/u, obtained from our sim- 
ulations, agrees with the order of magnitude with the 



5 



simple estimates by using the unipolar inductor model 
for the strongly magnetized star with a weakly magne- 
tized companion in a free space magnetosphere. This 
leads to similar effective temperatures, when the same 
assumptions (optically thick magnetosphere) are consid- 
ered, in the hard X-ray of the spectrum. Notice that in 
the extreme case where the primary star (for the U/u 
case) is a magnetar of B ^ 10^^ G the luminosity will 
increase to Lu/u ^ 1-6 x 10^^ erg/s, while that the black 
body temperature will rise up to Tu/^ = 2 x 10^ K, well 
inside the 7-ray range. 

We have shown global magnetic fields within a bi- 
nary system can give rise to a rich phenomenology that 
powers strong emissions on the electromagnetic side {o:^ 
10^^~^^(5/10"'^"'^G)^ erg/s) prior to the merger, setting 
ideal conditions for several emission processes. These lu- 
minosities are at the level of the brightest pulsars and 
beyond and would bear particular characteristics tied to 
the orbital behavior. This, together with the power emit- 
ted in gravitational waves indicate the system is strongly 
radiative in multiple bands and channels. We have also 
identified possible features that can lead to observable 
signals tied to the orbital behavior of the system. The 
consequential time-dependent nature of possible emis- 
sions might help in their detection, especially if some 
prior localization (in time and space) are provided by 
gravitational wave information. 
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